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ABSTRACT 

I present a simple scheme for the treatment of gravitational interactions on galactic scales. In analogy 
to known mechanisms of quantum field theory, I assume ad hoc that gravitation is mediated by virtual 
exchange particles - gravitons - with very small but non-zero masses. The resulting density and mass 
profiles are proportional to the mass of the gravitating body. The mass profile scales with the centripetal 
acceleration experienced by a test particle orbiting the central mass; this comes at the cost of postulating 
a universal characteristic acceleration oq ps 4.3 x 10~ 12 ms~ 2 (or 87rao « 1.1 x 10~ 10 ms -2 ). The 
scheme predicts the asymptotic flattening of galactic rotation curves, the Tully-Fisher/Faber- Jackson 
relations, the mass discrepancy-acceleration relation of galaxies, the surface brightness-acceleration 
relation of galaxies, the kinematics of galaxy clusters, and "Renzo's rule" correctly; additional (dark) 
mass components are not required. Given that it is based on various ad-hoc assumptions, and given 
further limitations, the scheme I present is not yet a consistent theory of gravitation; rather, it is a "toy 
model" providing a convenient scaling law that simplifies the description of gravity on galactic scales. 

Key words : Gravitation — Galaxies: kinematics and dynamics 



1. INTRODUCTION 

It is well-known that the treatment of gravitational in- 
teraction on galactic scales comes with a serious flaw 
denoted as the missing mass problem. Since the semi- 
nal works by Zwicky (1933) and Rubin et al. (1980), it 
has become clear that the dynamical mass necessary to 
explain the kinematics of galactic systems exceeds the 
actually observed luminous mass by up to one order 
of magnitude. On the one hand, this observation has 
lead to the postulate of non-luminous and non-baryonic 
dark matter that accounts for the observed discrep- 
ancy (Ostriker & Peebles 1973), eventually evolving 
into the modern ACDM standard model of cosmology 
(e.g., Bahcall et al. 1999). On the other hand, various 
proposals for modified laws of inertia and/or gravity 
have been made to address the missing mass problem, 
also in view of open questions regarding the standard 
cosmology (e.g. Kroupa 2012); the most successful can- 
didate to date appears to be "Milgrom's law" of modi- 
fied inertia (Milgrom 1983; Famaey & McGaugh 2012). 

To date, a fully satisfactory description of gravi- 
tational interaction on galactic scales has not been 
found. Dark matter models are struggling with "fine- 
tuning problems" regarding the distribution of dark 
matter within and around galaxies (e.g. Kroupa 2012), 
whereas "Milgrom's law" is based on an - entirely em- 
pirical - ad-hoc modification of Newton's law of inertia 
(e.g. Famaey & McGaugh 2012). 

My work aims at a simple, convenient, physically 



motivated description of gravity on galactic scales. In 
analogy to the well-established concepts of particle 
physics and quantum field theories (see e.g. Griffith 
2008 for an overview), I assume ad hoc that gravita- 
tional interactions between masses are mediated by dis- 
crete virtual particles dubbed gravitons. Assuming fur- 
ther that these gravitons have non-zero mass and obey 
certain rules of interaction, the resulting total mass dis- 
tributions (composed of source masses plus graviton 
distributions) can be treated by standard Newtonian 
dynamics (see e.g. Binney & Tremaine 2008 for a thor- 
ough review of stellar and galactic dynamics). The re- 
sulting scaling relations for mass profiles and rotation 
curves agree satisfactorily with observations. 

2. GRAVITATIONAL INTERACTION 
2.1 Fundamental Assumptions 

In order to provide a scheme of gravitational interac- 
tion, I start with the following fundamental assump- 
tions: 

Al. Gravitational interactions are mediated by dis- 
crete particles, gravitons. 

A2. Gravitons are virtual particles arising from quan- 
tum fluctuations - the mechanism commonly em- 
ployed in quantum field theories. 

A3. Gravitons have a non-zero mass. 

A4. Interactions may only occur (1) between two real 
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masses; (2) between a real mass and a graviton 
emitted by a real mass. 

2.2 Consequences 

Assumptions Al-4 lead to several consequences that, 
partially, can be tested by observations. Unless stated 
otherwise, the terms "gravitational system" or "dy- 
namical system" denote a system composed of a point- 
like, real, baryonic, luminous source mass Mq and a 
quasi-massless test particle orbiting M on a circular 
orbit at radial distance r with circular speed v c . 

2.2.1 Gravitons Are "Dark" 

Postulating gravitons as mediators of gravitational in- 
teractions (assumption Al) implies that they couple to 
mass, with mass playing the role of the "charge" of 
gravity. Gravitons should not be expected to take part 
in electromagnetic/electroweak or strong nuclear inter- 
actions (sec e.g. Griffith 2008 for an overview). Ac- 
cordingly, gravitons do not interact with photons elec- 
tromagnetically, they are "dark" in the sense of being 
"invisible" to electromagnetic radiation. 

2.2.2 Graviton Masses 

The creation of virtual massive particles from quantum 
fluctuations (assumptions A2 and A3) is, a priori, a vi- 
olation of the principle of conservation of energy. This 
violation is only possible as long as it is "temporary" 
within the limits of Heisenberg's uncertainty relation 
between uncertainties in energy and time, meaning here 
specifically 

AEAt^Ti (1) 

where AE = m g c 2 is the energy of a graviton with 
mass m g , c is the speed of light, At is the lifetime of 
the virtual particle, and Ti — 1.055 x 10~ 34 Js is the 
reduced form of Planck's constant (e.g. Griffith 2008). 

A range of gravitational interaction on cosmological 
scales requires graviton lifetimes on the order of the 
Hubble time, meaning At J> 1.4xl0 10 yrs. Accordingly, 
we require 



m 9 = — £ 3 x !0~ 69 kg« 1.5 x 10- 33 cVc- 2 . (2) 

This value is about 40 orders of magnitude smaller than 
the masses of virtual exchange particles occurring in 
particle physics - a necessary consequence of the vast 
difference in size scales (see also Goldhaber & Nicto 
2010 for a recent review of graviton mass limits). 

2.2.3 Density and Mass Profiles 

I assume in the following that the source mass M radi- 
ates away gravitons without loosing mass itself; this is a 



corollary of assumption A2. Consistency with the clas- 
sical theory of gravitational fields requires that (1) the 
graviton density profile scales linearly with Mo, and (2) 
the density profile obeys the inverse-square-of-distance 
law of gravitational fields (and of radiation fields in 
general). Accordingly, we find 

p = M Q (3r- 2 (3) 

where p is the mass density of the graviton distribution, 
r is the radial distance from Mo, and (3 is a scaling 
parameter of the dimension of an inverse length. 

The mass profile is found from integration of p(r). 
Demanding that /? does not depend on r explicitly,* we 
find the extra mass enclosed within a distance r to be 

M cx = 4TrM f3r (4) 
whereas the total enclosed mass is given by 

M tot - M + M cx = M (1 + 4tt r) . (5) 

2.2.4 Enclosed Mass vs. Centripetal Acceleration 

To permit a comparison of arbitrary gravitational sys- 
tems, we need to bring the scaling parameter ft in Eq. 4 
into a form that (1) provides a characteristic scale for 
a given dynamical system, and (2) does not depend on 
the coordinate r explicitly. Both conditions are fulfilled 
by choosing 



where v c 2 /2 is the kinetic energy per unit mass of a 
test particle orbiting M with circular speed v c , and ao 
is a constant of the dimension of an acceleration. Con- 
veniently, the kinetic energy is scale-free, i.e. invariant 
under simultaneous scalings of coordinates and time by 
factors x 7^ 0. Inserting /? into Eqs. 4 and 5 leads to a 
new expression for the enclosed total mass, 

M tot = M + 8tt^ (7) 

with a c = v c 2 /r being the centripetal acceleration ex- 
perienced by a test particle orbiting M with circular 
speed v c at radial distance r . 

2.2.5 Absence of Self-Interaction 

Assumption A4 is crucial for the self-consistency of the 
"graviton picture" of gravitation I use. Under assump- 
tions Al and A3, gravity is mediated by gravitons that 
have mass themselves. As mass can be regarded as the 



* Meaning here specifically: f} may be a function of a given well- 
defined distance r' but not of the coordinate r in general. This 
permits integration of p(r) over r from to r' and a subsequent 
substitution r' — > r. 
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"charge" of gravitation, we may, a priori, expect in- 
teractions also between virtual masses, leading to the 
creation of additional masses, that interact again, et 
cetera. Those interactions of higher order can result in 
a divergence in the creation of mediator particles, lead- 
ing to massive modifications of the overall interactions 
- especially to violations of the r~ 2 law of gravity (and 
radiation in general). This is a situation commonly en- 
countered - and treated - in quantum chromodynam- 
ics, expressed especially in the confinement of quarks 
(Griffith 2008). 

Assuming a quasi-massless test particle interacting 
with Mo, assumption A4 prevents those divergences: 

• By permitting interactions between real masses, 
we ensure that the test mass "sees" the source 
mass Mo in the usual - Newtonian - way. 

• By permitting interactions between real masses 
and virtual masses emitted by a real mass, we 
ensure that the test mass "sees" the distribution 
of primary gravitons emitted by Mo- 

• By preventing any other interaction, we suppress 
interactions of higher order, potentially leading to 
divergences and/or confinements. 

As mentioned above, the conditions imposed by as- 
sumption A4 are far from obvious. In quantum field 
theories, those conditions are usually traced back to 
conserved quantities that impose selection rules for al- 
lowed interactions; prominent examples are spins or 
the color charges of quantum chromodynamics (Grif- 
fith 2008). 

3. TESTS 

3.1 Galactic Rotation Curves 

The circular speed of a test particle orbiting M is re- 
lated to Eq. 7 like 

^ = ^ = ^(1 + 8^ (8) 
r r V a cJ 

with G denoting Newton's constant. As long as 
87rao/a c <C 1, the first summand in the bracket domi- 
nates, resulting in (quasi-)Keplerian motion. However, 
in the limit 8nao/a c ^> 1, the second summand domi- 
nates and we have 



8tt G M a Q 
r a c 



(9) 



Exploiting the fact that a c = v c 2 /r, we can rewrite 
Eq. 9 like 



w c 4 « 87r G M ao = const. 



(10) 



Accordingly, one finds that, in the limit of strong cen- 
tripetal accelerations, the circular speed is Keplerian; 



it approaches a constant value when approaching the 
limit of weak centripetal accelerations. This is in agree- 
ment with the observed behavior of galactic rotation 
curves (Rubin et al. 1980). 

3.2 Tully-Fisher/Faber-Jackson Relations 

Next to the fact that rotation curves flatten out asymp- 
totically, Eq. 10 provides additional information: it 
states that, asymptotically, 



v c 4 oc M 



(11) 



This is in agreement with the empirical Tully-Fisher 
relation between rotation speeds of spiral galaxies and 
their luminous masses (Tully & Fisher 1977). A re- 
lation as predicted by Eq. 10 is observationally estab- 
lished over ten orders of magnitude of mass for gravita- 
tionally bound systems, from dwarf galaxies to galaxy 
clusters (cf. Fig. 48 of Famaey & McGaugh 2012). 
Likewise, Eq. 10 provides an explanation for the em- 
pirical Faber-, Jackson relation between the luminosity 
- which is proportional to the source mass in my inter- 
action scheme - of elliptical galaxies and their velocity 
dispersions (Faber & Jackson 1976). 

3.3 Mass Discrepancy Acceleration Relation 

Equation 7 links the total enclosed mass, the source 
mass, and the centripetal acceleration experienced by a 
test particle orbiting the source mass. By rewriting this 
equation slightly, we can formulate a mass discrepancy- 
acceleration relation (McGaugh 2004) that links the 
mass discrepancy, defined as the ratio M to t/M , and 
acceleration directly: 



^ = 1 + 8--- 

M a c 



(12) 



Relation 12 is illustrated in Fig. la for a range of 
weak centripetal accelerations a c . In addition, Fig. lb 
shows the mass discrepancy as function of the New- 
tonian gravitational acceleration expected if Mo were 
equal to the total enclosed mass, = a c /(M to t/ Mq). 
Comparison to Fig. 10 of Famaey & McGaugh (2012) 
indicates a satisfactory agreement between Eq. 12 and 
observations for a = 4.3 x 10~ 12 ms~ 2 (±5%). 

3.4 Surface Brightness vs. Acceleration 

In the limit of weak centripetal accelerations, we can 
use Eq. 10 and rewrite it like 



8nGao 



Mo 



(13) 



For a sufficiently flat stellar system, we can define a 
(luminous) surface density £ oc Mo/r 2 . Inserting E 
into Eq. 13 and exploiting the fact that a c — v c 2 /r, we 
may conclude that 
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Fig. 1. — Theoretical mass discrepancy-acceleration re- 
lations, a. The ratio M to t/Mo as function of centripetal 
acceleration a c according to Eq. 12. b. M to t/Mo as func- 
tion of the Newtonian gravitational acceleration expected 
if Mo were the total mass, gw = a c /(M to t/Mo). The char- 
acteristic acceleration used here is ao = 4.3 x 10~ 12 ms~ 2 . 
These diagrams should be compared to Fig. 10 of Famaey 
& McGaugh (2012). 



a c 2 oc £ or a c oc S 1/2 . (14) 

Such a correlation is indeed observed (see Fig. 9 of 
Famaey & McGaugh 2012). 

3.5 Galaxy Clusters 

The interaction scheme I use should not be restricted 
to galaxies but also become visible in the dynamics of 
galaxy clusters, likewise leading to observable mass dis- 
crepancies. Unfortunately, application of Eq. 7 is not 
straightforward because galaxy clusters are not point- 
like masses but extended mass distributions with most 
of their masses apparently stored in the X-ray luminous 



intra-cluster medium (e.g. Sparke & Gallagher 2007). 

I estimate the mass discrepancy crudely - accurate 
to factors of a few - by using the velocity dispersions 
and core radii of the clusters as proxies for v c and 
orbit radii, respectively. For the four clusters Virgo, 
Coma, Perseus, and RDCS 1252.9-2927 (from Table 
7.1 of Sparke & Gallagher 2007), this leads to theoret- 
ical mass discrepancies of about 2 to 4; the actually 
observed values are about 6 to 10. Given the crudc- 
ness of the estimate, I conclude that theoretical and 
observational values are consistent with each other; es- 
pecially, both are of the same order of magnitude <;10. 

3.6 "Renzo's Rule" 

Observationally, dynamical and (luminous) baryonic 
surface densities in disk galaxies are closely related, 
expressed - among others - by "Renzo's rule" (Sancisi 
2004) which states that "for any feature in the lumi- 
nosity profile, there is a corresponding feature in the 
rotation curve" . This empirical rule has been used as 
an argument against the presence of dark matter dis- 
tributions in disk galaxies because these are, a priori, 
independent from structures in the distributions of lu- 
minous mass (Famaey & McGaugh 2012). 

In the frame of the interaction scheme presented in 
this work, "Renzo's rule" is self-evident due to the 
unique relation between (luminous) source mass and 
total enclosed mass (Eqs. 7, 8). 

4. DISCUSSION 

Starting from the ad-hoc assumptions Al-4, I arrive 
at a simple description of gravitational interaction on 
galactic scales. This "graviton picture" of gravity 
makes several predictions about galactic dynamics that 
are consistent with observations. Compared to stan- 
dard models of galactic dynamics, the scheme I pro- 
pose seems preferable because it removes the need for 
a separate "dark mass" component - thus providing 
a simplified description. Notably, the "graviton pic- 
ture" seems to be the only model to date that predicts 
the mass discrepancy-acceleration relation (Sect. 3.3, 
Fig. 1), and thus the transitional regime between New- 
tonian and modified gravitational interaction, correctly 
- modified laws of gravity/inertia or dark matter mod- 
els in general make predictions on the asymptotic be- 
havior of stellar systems only. 

Although it seems that my work is the first to ex- 
plicitly apply the assumption of massive gravitons to 
galactic dynamics, the concept of massive gravitons - 
or "massive gravity" - itself has been discussed in field 
theories for more than seven decades (see, e.g., Gold- 
haber & Nieto 2010; Hinterbichler 2012 for reviews) 
and has been applied to various aspects of cosmology 
(e.g. Cardone, Radicella & Parisi 2012). 

Even though my interaction scheme appears success- 
ful on galactic scales, it is obviously incomplete. It is 



GRAVITATION ON GALACTIC SCALES 



5 



formulated in the frame of classical Newtonian dynam- 
ics, meaning it is intrinsically non-relativistic. This 
implies that the scheme breaks down at small scales 
where relativistic deviations from Newtonian dynamics 
become noticeable - on scales of the size of the solar 
system and less - as well as at large scales where cos- 
mological effects, like the Hubble flow, become evident 
on scales of mega-parsecs and beyond. 

A further peculiarity of the description I provide 
is encoded in the formulation of the scaling param- 
eter /? in Eq. 6.^ The choice I make is, to a cer- 
tain extent, technically self-evident: Eq. 6 provides 
a characteristic kinematic size scale, is independent 
of the radial coordinate, and is based on the scale- 
invariant, conserved kinetic energy. However, the ap- 
parent astrophysical consequences discussed through- 
out Sect. 3 enforce a remarkable postulate: there is a 
characteristic acceleration ao ~ 4.3 x 10~ 12 ms~ 2 (or 
87rao ~ 1.1 x 10~ 10 ms~ 2 ) which is a universal constant 
of nature. The origin of this characteristic acceleration 
remains unclear;-' 1 I note however that it is well-known 
empirically (cf. e.g. Famaey & McGaugh 2012, and 
references therein). 

Given that the interaction scheme I propose is based 
on various ad-hoc assumptions, and given its limita- 
tions discussed above, it is clear that the "graviton pic- 
ture" is not yet a consistent theory of gravitation. In- 
stead, it should be seen as a "toy model" that provides 
a simple scaling law for gravitational interaction on 
galactic scales and that is simpler than standard models 
(that require additional mass components) . This seems 
valuable especially in view of modern numerical simu- 
lations exploring galactic to cosmological scales (e.g. 
Kim et al. 2011). 

Despite the various limitations stated above, one can 
think of several observational signatures that could be 
used to further probe the "graviton picture" , in ad- 
dition to the "usual suspects" discussed throughout 
Sect. 3: 

Dynamical friction. To a large extent, the dynamics 
of stellar systems is governed by dynamical friction (cf. 
Ch. 8.1 of Binney & Tremaine 2008). The assumption 
of dark matter halos around galaxies predicts a variety 
of dynamical signatures ranging from the evolution of 
binary star orbits (e.g. Hernandez & Lee 2008) to the 
effects of galaxy encounters (e.g. Dubinski, Mihos & 
Hcrnquist 1999); accordingly, models for dark matter 
and/or modified laws of gravity can, and have been, 
examined by comparing the results of N-body simula- 



tl also note an analogy provided by atomic physics. Rydbcrg's 
constant Roc, is the only constant of nature of the dimension of an 
inverse length (e.g. Lang 2006). In Bohr's model, Roc is the in- 
verse wavelength of a photon with the energy of the ground state 
of a hydrogen atom. This means that Roc, actually provides an 
energy scale; it is only a constant of nature because all hydrogen 
atoms arc identical. Likewise, fi provides a characteristic energy 
scale for dynamical systems. 

* Though this discussion might be loosely related to the discussion 
of Mach's principle; see e.g. Bondi & Samuel (1997). 



tions with observations of those systems. The "gravi- 
ton picture" makes a strong prediction: as a corollary 
of assumption A4, only baryonic masses experience dy- 
namical friction in a "graviton halo" ; the graviton dis- 
tributions themselves do not interact. This is sharply 
distinct from the case of dark matter halos where both 
baryonic and dark mass particles are supposed to ex- 
perience dynamical friction. Likewise, modified laws 
of gravity do not introduce any extra masses, mean- 
ing they do not predict dynamical friction beyond the 
one due to luminous mass. Accordingly, observational 
evidence, e.g. in galaxy collisions, for interactions cor- 
responding to dynamical friction of baryonic mass in 
"dark halos" without interaction among those halos 
would provide support for the "graviton picture" . 
Collisions of galaxy clusters. The "graviton pic- 
ture" demands a direct coupling of a luminous bary- 
onic source mass and the extra mass due to the sur- 
rounding graviton distribution (cf. Eqs. 4, 5, 7). In 
the standard dark matter picture, luminous and dark 
mass may be distributed independently A test of these 
scenarios is provided by colliding clusters of galaxies: 
when the cluster cores pass through each other, the hot, 
X-ray bright intra-cluster gas, which comprises a large 
fraction of the baryonic cluster masses, experiences col- 
lisional ram-pressure and lags behind the collisionless 
galaxies and dark matter halos - if those actually ex- 
ist. Accordingly, a comparison of the spatial distribu- 
tions of total mass - derived from gravitational lensing 
of background sources - and luminous mass is able to 
test dark matter and/or modified gravity models. An 
important example is provided by the "Train- Wreck 
Cluster" A520 (Markevitch et al. 2005). In A520, ob- 
servations find that the spatial distributions of total 
and luminous mass coincide (Jee et al. 2012) - in agree- 
ment with the "graviton picture" . Another example is 
provided by the "Bullet Cluster" 1E0657-56 (Barrena 
et al. 2002). Here the situation is more complicated: 
on the one hand, observations indicate a spatial sepa- 
ration of dark and luminous mass; on the other hand, 
the cluster kinematics is found to be inconsistent with 
ACDM cosmology (Lee & Komatsu 2010). 
Gravitational lensing. The "graviton picture" makes 
a specific prediction for the distribution of the non- 
luminous extra mass around a luminous source mass: 
as it assumes that gravity is the result of a radiation 
composed of massive particles, the density distribution 
of these particles has to obey the usual r -2 law of ra- 
diation (cf. Eq. 3). This prediction has been tested by 
gravitational lensing studies that probe the density pro- 
files of galaxies; observational results are in agreement 
with p oc r~ 2 as demanded by the "graviton picture" 
(Gavazzi et al. 2007). 

Cosmology. A non-zero mass implies a limited life 
time of gravitons (see Sect. 2.2.2). Accordingly, a decay 
of gravitons on cosmological time scales should modify 
the cosmic expansion history; especially, an accelera- 
tion corresponding to the one commonly ascribed to 
dark energy (Bahcall et al. 1999) might occur. Recent 
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studies (Cardone, Radicella & Parisi 2012) indeed con- 
clude that massive gravity is able to reproduce the ob- 
served cosmological parameters. 

Evidently, further theoretical as well as observa- 
tional work is required to understand the limitations 
and implications of the interaction scheme I present in 
this study. Eventually the "graviton picture" might 
play a role for galactic dynamics which is analogous to 
the one played by Bohr's "planetary system" model of 
the hydrogen atom for atomic physics: even though it is 
an over-simplified, coarse approximation of the actual 
physical situation, it provides a successful and conve- 
nient quantitative description for a wide range of phe- 
nomena. Ironically, the "graviton picture" resolves the 
ambiguity "dark matter or modified gravity" by assum- 
ing "dark matter through modified gravity" . 

5. CONCLUSIONS 

Based on a set of four ad-hoc assumptions, I present a 
simple scheme for the treatment of gravitational inter- 
action on galactic scales: 

1. I assume that gravitation is mediated by gravi- 
tons, which are virtual, discrete exchange particles 
with non-zero mass. These gravitons ought to in- 
teract gravitationally but not electromagnetically, 
meaning they are "invisible" to light. 

2. Demanding consistency with classical field theory, 
the mass and density profiles of the graviton dis- 
tribution around a source mass are proportional 
to the mass of the source. The density profile 
obeys the usual inverse-square-of-distance law of 
gravitational fields. 

3. The mass profile contains a scaling factor f3 that 
can be written as the ratio of a universal, con- 
stant acceleration and the kinetic energy per unit 
mass of a test particle orbiting the source mass. 
Empirically, the characteristic acceleration is a ~ 
4.3 x 10~ 12 ms~ 2 . In the chosen formulation, the 
mass profile is a function of the centripetal accel- 
eration experienced by a test particle orbiting the 
source mass. 

4. Mass profiles and resulting rotation curves pre- 
dict the asymptotic flattening of galactic rota- 
tion curves, the Tully-Fisher/Faber- Jackson rela- 
tions, the mass discrepancy-acceleration relation 
of galaxies, the surface brightness-acceleration re- 
lation of galaxies, the kinematics of galaxy clus- 
ters, and "Renzo's rule" correctly. Additional fu- 
ture observational tests could be provided by dy- 
namical friction in stellar systems, collisions of 
galaxy clusters, gravitational lensing by galaxies, 
and the evolution of cosmic expansion. 

5. The interaction scheme is non-relativistic, mean- 
ing it should break down on small - solar system 
- as well as cosmological scales. 



Within its known limitations, the "graviton picture" 
provides a useful tool for studies of galactic dynamics. 
Obviously, its range of validity needs to be explored 
more carefully in future observational as well as theo- 
retical/numerical studies. 
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